3.271 \(\int \frac{x^4 (c+d x^2)^2}{(a+b x^2)^2} \, dx\)

Optimal. Leaf size=145 \[ \frac{x^5 (b c-a d)^2}{2 a b^2 \left (a+b x^2\right )}-\frac{x^3 (3 b c-7 a d) (b c-a d)}{6 a b^3}+\frac{x (3 b c-7 a d) (b c-a d)}{2 b^4}-\frac{\sqrt{a} (3 b c-7 a d) (b c-a d) \tan ^{-1}\left (\frac{\sqrt{b} x}{\sqrt{a}}\right )}{2 b^{9/2}}+\frac{d^2 x^5}{5 b^2} \]

[Out]

((3*b*c - 7*a*d)*(b*c - a*d)*x)/(2*b^4) - ((3*b*c - 7*a*d)*(b*c - a*d)*x^3)/(6*a*b^3) + (d^2*x^5)/(5*b^2) + ((
b*c - a*d)^2*x^5)/(2*a*b^2*(a + b*x^2)) - (Sqrt[a]*(3*b*c - 7*a*d)*(b*c - a*d)*ArcTan[(Sqrt[b]*x)/Sqrt[a]])/(2
*b^(9/2))

________________________________________________________________________________________

Rubi [A]  time = 0.13212, antiderivative size = 145, normalized size of antiderivative = 1., number of steps used = 5, number of rules used = 4, integrand size = 22, \(\frac{\text{number of rules}}{\text{integrand size}}\) = 0.182, Rules used = {463, 459, 302, 205} \[ \frac{x^5 (b c-a d)^2}{2 a b^2 \left (a+b x^2\right )}-\frac{x^3 (3 b c-7 a d) (b c-a d)}{6 a b^3}+\frac{x (3 b c-7 a d) (b c-a d)}{2 b^4}-\frac{\sqrt{a} (3 b c-7 a d) (b c-a d) \tan ^{-1}\left (\frac{\sqrt{b} x}{\sqrt{a}}\right )}{2 b^{9/2}}+\frac{d^2 x^5}{5 b^2} \]

Antiderivative was successfully verified.

[In]

Int[(x^4*(c + d*x^2)^2)/(a + b*x^2)^2,x]

[Out]

((3*b*c - 7*a*d)*(b*c - a*d)*x)/(2*b^4) - ((3*b*c - 7*a*d)*(b*c - a*d)*x^3)/(6*a*b^3) + (d^2*x^5)/(5*b^2) + ((
b*c - a*d)^2*x^5)/(2*a*b^2*(a + b*x^2)) - (Sqrt[a]*(3*b*c - 7*a*d)*(b*c - a*d)*ArcTan[(Sqrt[b]*x)/Sqrt[a]])/(2
*b^(9/2))

Rule 463

Int[((e_.)*(x_))^(m_.)*((a_) + (b_.)*(x_)^(n_))^(p_)*((c_) + (d_.)*(x_)^(n_))^2, x_Symbol] :> -Simp[((b*c - a*
d)^2*(e*x)^(m + 1)*(a + b*x^n)^(p + 1))/(a*b^2*e*n*(p + 1)), x] + Dist[1/(a*b^2*n*(p + 1)), Int[(e*x)^m*(a + b
*x^n)^(p + 1)*Simp[(b*c - a*d)^2*(m + 1) + b^2*c^2*n*(p + 1) + a*b*d^2*n*(p + 1)*x^n, x], x], x] /; FreeQ[{a,
b, c, d, e, m, n}, x] && NeQ[b*c - a*d, 0] && IGtQ[n, 0] && LtQ[p, -1]

Rule 459

Int[((e_.)*(x_))^(m_.)*((a_) + (b_.)*(x_)^(n_))^(p_.)*((c_) + (d_.)*(x_)^(n_)), x_Symbol] :> Simp[(d*(e*x)^(m
+ 1)*(a + b*x^n)^(p + 1))/(b*e*(m + n*(p + 1) + 1)), x] - Dist[(a*d*(m + 1) - b*c*(m + n*(p + 1) + 1))/(b*(m +
 n*(p + 1) + 1)), Int[(e*x)^m*(a + b*x^n)^p, x], x] /; FreeQ[{a, b, c, d, e, m, n, p}, x] && NeQ[b*c - a*d, 0]
 && NeQ[m + n*(p + 1) + 1, 0]

Rule 302

Int[(x_)^(m_)/((a_) + (b_.)*(x_)^(n_)), x_Symbol] :> Int[PolynomialDivide[x^m, a + b*x^n, x], x] /; FreeQ[{a,
b}, x] && IGtQ[m, 0] && IGtQ[n, 0] && GtQ[m, 2*n - 1]

Rule 205

Int[((a_) + (b_.)*(x_)^2)^(-1), x_Symbol] :> Simp[(Rt[a/b, 2]*ArcTan[x/Rt[a/b, 2]])/a, x] /; FreeQ[{a, b}, x]
&& PosQ[a/b]

Rubi steps

\begin{align*} \int \frac{x^4 \left (c+d x^2\right )^2}{\left (a+b x^2\right )^2} \, dx &=\frac{(b c-a d)^2 x^5}{2 a b^2 \left (a+b x^2\right )}-\frac{\int \frac{x^4 \left (-2 b^2 c^2+5 (b c-a d)^2-2 a b d^2 x^2\right )}{a+b x^2} \, dx}{2 a b^2}\\ &=\frac{d^2 x^5}{5 b^2}+\frac{(b c-a d)^2 x^5}{2 a b^2 \left (a+b x^2\right )}-\frac{((3 b c-7 a d) (b c-a d)) \int \frac{x^4}{a+b x^2} \, dx}{2 a b^2}\\ &=\frac{d^2 x^5}{5 b^2}+\frac{(b c-a d)^2 x^5}{2 a b^2 \left (a+b x^2\right )}-\frac{((3 b c-7 a d) (b c-a d)) \int \left (-\frac{a}{b^2}+\frac{x^2}{b}+\frac{a^2}{b^2 \left (a+b x^2\right )}\right ) \, dx}{2 a b^2}\\ &=\frac{(3 b c-7 a d) (b c-a d) x}{2 b^4}-\frac{(3 b c-7 a d) (b c-a d) x^3}{6 a b^3}+\frac{d^2 x^5}{5 b^2}+\frac{(b c-a d)^2 x^5}{2 a b^2 \left (a+b x^2\right )}-\frac{(a (3 b c-7 a d) (b c-a d)) \int \frac{1}{a+b x^2} \, dx}{2 b^4}\\ &=\frac{(3 b c-7 a d) (b c-a d) x}{2 b^4}-\frac{(3 b c-7 a d) (b c-a d) x^3}{6 a b^3}+\frac{d^2 x^5}{5 b^2}+\frac{(b c-a d)^2 x^5}{2 a b^2 \left (a+b x^2\right )}-\frac{\sqrt{a} (3 b c-7 a d) (b c-a d) \tan ^{-1}\left (\frac{\sqrt{b} x}{\sqrt{a}}\right )}{2 b^{9/2}}\\ \end{align*}

Mathematica [A]  time = 0.0892508, size = 138, normalized size = 0.95 \[ \frac{x \left (3 a^2 d^2-4 a b c d+b^2 c^2\right )}{b^4}-\frac{\sqrt{a} \left (7 a^2 d^2-10 a b c d+3 b^2 c^2\right ) \tan ^{-1}\left (\frac{\sqrt{b} x}{\sqrt{a}}\right )}{2 b^{9/2}}+\frac{2 d x^3 (b c-a d)}{3 b^3}+\frac{a x (b c-a d)^2}{2 b^4 \left (a+b x^2\right )}+\frac{d^2 x^5}{5 b^2} \]

Antiderivative was successfully verified.

[In]

Integrate[(x^4*(c + d*x^2)^2)/(a + b*x^2)^2,x]

[Out]

((b^2*c^2 - 4*a*b*c*d + 3*a^2*d^2)*x)/b^4 + (2*d*(b*c - a*d)*x^3)/(3*b^3) + (d^2*x^5)/(5*b^2) + (a*(b*c - a*d)
^2*x)/(2*b^4*(a + b*x^2)) - (Sqrt[a]*(3*b^2*c^2 - 10*a*b*c*d + 7*a^2*d^2)*ArcTan[(Sqrt[b]*x)/Sqrt[a]])/(2*b^(9
/2))

________________________________________________________________________________________

Maple [A]  time = 0.009, size = 196, normalized size = 1.4 \begin{align*}{\frac{{d}^{2}{x}^{5}}{5\,{b}^{2}}}-{\frac{2\,{x}^{3}a{d}^{2}}{3\,{b}^{3}}}+{\frac{2\,c{x}^{3}d}{3\,{b}^{2}}}+3\,{\frac{{a}^{2}{d}^{2}x}{{b}^{4}}}-4\,{\frac{acdx}{{b}^{3}}}+{\frac{{c}^{2}x}{{b}^{2}}}+{\frac{{a}^{3}x{d}^{2}}{2\,{b}^{4} \left ( b{x}^{2}+a \right ) }}-{\frac{{a}^{2}cdx}{{b}^{3} \left ( b{x}^{2}+a \right ) }}+{\frac{ax{c}^{2}}{2\,{b}^{2} \left ( b{x}^{2}+a \right ) }}-{\frac{7\,{a}^{3}{d}^{2}}{2\,{b}^{4}}\arctan \left ({bx{\frac{1}{\sqrt{ab}}}} \right ){\frac{1}{\sqrt{ab}}}}+5\,{\frac{{a}^{2}cd}{{b}^{3}\sqrt{ab}}\arctan \left ({\frac{bx}{\sqrt{ab}}} \right ) }-{\frac{3\,a{c}^{2}}{2\,{b}^{2}}\arctan \left ({bx{\frac{1}{\sqrt{ab}}}} \right ){\frac{1}{\sqrt{ab}}}} \end{align*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

int(x^4*(d*x^2+c)^2/(b*x^2+a)^2,x)

[Out]

1/5*d^2*x^5/b^2-2/3/b^3*x^3*a*d^2+2/3/b^2*x^3*c*d+3/b^4*a^2*d^2*x-4/b^3*c*a*d*x+1/b^2*c^2*x+1/2*a^3/b^4*x/(b*x
^2+a)*d^2-a^2/b^3*x/(b*x^2+a)*c*d+1/2*a/b^2*x/(b*x^2+a)*c^2-7/2*a^3/b^4/(a*b)^(1/2)*arctan(b*x/(a*b)^(1/2))*d^
2+5*a^2/b^3/(a*b)^(1/2)*arctan(b*x/(a*b)^(1/2))*c*d-3/2*a/b^2/(a*b)^(1/2)*arctan(b*x/(a*b)^(1/2))*c^2

________________________________________________________________________________________

Maxima [F(-2)]  time = 0., size = 0, normalized size = 0. \begin{align*} \text{Exception raised: ValueError} \end{align*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate(x^4*(d*x^2+c)^2/(b*x^2+a)^2,x, algorithm="maxima")

[Out]

Exception raised: ValueError

________________________________________________________________________________________

Fricas [A]  time = 1.64322, size = 851, normalized size = 5.87 \begin{align*} \left [\frac{12 \, b^{3} d^{2} x^{7} + 4 \,{\left (10 \, b^{3} c d - 7 \, a b^{2} d^{2}\right )} x^{5} + 20 \,{\left (3 \, b^{3} c^{2} - 10 \, a b^{2} c d + 7 \, a^{2} b d^{2}\right )} x^{3} + 15 \,{\left (3 \, a b^{2} c^{2} - 10 \, a^{2} b c d + 7 \, a^{3} d^{2} +{\left (3 \, b^{3} c^{2} - 10 \, a b^{2} c d + 7 \, a^{2} b d^{2}\right )} x^{2}\right )} \sqrt{-\frac{a}{b}} \log \left (\frac{b x^{2} - 2 \, b x \sqrt{-\frac{a}{b}} - a}{b x^{2} + a}\right ) + 30 \,{\left (3 \, a b^{2} c^{2} - 10 \, a^{2} b c d + 7 \, a^{3} d^{2}\right )} x}{60 \,{\left (b^{5} x^{2} + a b^{4}\right )}}, \frac{6 \, b^{3} d^{2} x^{7} + 2 \,{\left (10 \, b^{3} c d - 7 \, a b^{2} d^{2}\right )} x^{5} + 10 \,{\left (3 \, b^{3} c^{2} - 10 \, a b^{2} c d + 7 \, a^{2} b d^{2}\right )} x^{3} - 15 \,{\left (3 \, a b^{2} c^{2} - 10 \, a^{2} b c d + 7 \, a^{3} d^{2} +{\left (3 \, b^{3} c^{2} - 10 \, a b^{2} c d + 7 \, a^{2} b d^{2}\right )} x^{2}\right )} \sqrt{\frac{a}{b}} \arctan \left (\frac{b x \sqrt{\frac{a}{b}}}{a}\right ) + 15 \,{\left (3 \, a b^{2} c^{2} - 10 \, a^{2} b c d + 7 \, a^{3} d^{2}\right )} x}{30 \,{\left (b^{5} x^{2} + a b^{4}\right )}}\right ] \end{align*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate(x^4*(d*x^2+c)^2/(b*x^2+a)^2,x, algorithm="fricas")

[Out]

[1/60*(12*b^3*d^2*x^7 + 4*(10*b^3*c*d - 7*a*b^2*d^2)*x^5 + 20*(3*b^3*c^2 - 10*a*b^2*c*d + 7*a^2*b*d^2)*x^3 + 1
5*(3*a*b^2*c^2 - 10*a^2*b*c*d + 7*a^3*d^2 + (3*b^3*c^2 - 10*a*b^2*c*d + 7*a^2*b*d^2)*x^2)*sqrt(-a/b)*log((b*x^
2 - 2*b*x*sqrt(-a/b) - a)/(b*x^2 + a)) + 30*(3*a*b^2*c^2 - 10*a^2*b*c*d + 7*a^3*d^2)*x)/(b^5*x^2 + a*b^4), 1/3
0*(6*b^3*d^2*x^7 + 2*(10*b^3*c*d - 7*a*b^2*d^2)*x^5 + 10*(3*b^3*c^2 - 10*a*b^2*c*d + 7*a^2*b*d^2)*x^3 - 15*(3*
a*b^2*c^2 - 10*a^2*b*c*d + 7*a^3*d^2 + (3*b^3*c^2 - 10*a*b^2*c*d + 7*a^2*b*d^2)*x^2)*sqrt(a/b)*arctan(b*x*sqrt
(a/b)/a) + 15*(3*a*b^2*c^2 - 10*a^2*b*c*d + 7*a^3*d^2)*x)/(b^5*x^2 + a*b^4)]

________________________________________________________________________________________

Sympy [B]  time = 1.12611, size = 280, normalized size = 1.93 \begin{align*} \frac{x \left (a^{3} d^{2} - 2 a^{2} b c d + a b^{2} c^{2}\right )}{2 a b^{4} + 2 b^{5} x^{2}} + \frac{\sqrt{- \frac{a}{b^{9}}} \left (a d - b c\right ) \left (7 a d - 3 b c\right ) \log{\left (- \frac{b^{4} \sqrt{- \frac{a}{b^{9}}} \left (a d - b c\right ) \left (7 a d - 3 b c\right )}{7 a^{2} d^{2} - 10 a b c d + 3 b^{2} c^{2}} + x \right )}}{4} - \frac{\sqrt{- \frac{a}{b^{9}}} \left (a d - b c\right ) \left (7 a d - 3 b c\right ) \log{\left (\frac{b^{4} \sqrt{- \frac{a}{b^{9}}} \left (a d - b c\right ) \left (7 a d - 3 b c\right )}{7 a^{2} d^{2} - 10 a b c d + 3 b^{2} c^{2}} + x \right )}}{4} + \frac{d^{2} x^{5}}{5 b^{2}} - \frac{x^{3} \left (2 a d^{2} - 2 b c d\right )}{3 b^{3}} + \frac{x \left (3 a^{2} d^{2} - 4 a b c d + b^{2} c^{2}\right )}{b^{4}} \end{align*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate(x**4*(d*x**2+c)**2/(b*x**2+a)**2,x)

[Out]

x*(a**3*d**2 - 2*a**2*b*c*d + a*b**2*c**2)/(2*a*b**4 + 2*b**5*x**2) + sqrt(-a/b**9)*(a*d - b*c)*(7*a*d - 3*b*c
)*log(-b**4*sqrt(-a/b**9)*(a*d - b*c)*(7*a*d - 3*b*c)/(7*a**2*d**2 - 10*a*b*c*d + 3*b**2*c**2) + x)/4 - sqrt(-
a/b**9)*(a*d - b*c)*(7*a*d - 3*b*c)*log(b**4*sqrt(-a/b**9)*(a*d - b*c)*(7*a*d - 3*b*c)/(7*a**2*d**2 - 10*a*b*c
*d + 3*b**2*c**2) + x)/4 + d**2*x**5/(5*b**2) - x**3*(2*a*d**2 - 2*b*c*d)/(3*b**3) + x*(3*a**2*d**2 - 4*a*b*c*
d + b**2*c**2)/b**4

________________________________________________________________________________________

Giac [A]  time = 1.1463, size = 211, normalized size = 1.46 \begin{align*} -\frac{{\left (3 \, a b^{2} c^{2} - 10 \, a^{2} b c d + 7 \, a^{3} d^{2}\right )} \arctan \left (\frac{b x}{\sqrt{a b}}\right )}{2 \, \sqrt{a b} b^{4}} + \frac{a b^{2} c^{2} x - 2 \, a^{2} b c d x + a^{3} d^{2} x}{2 \,{\left (b x^{2} + a\right )} b^{4}} + \frac{3 \, b^{8} d^{2} x^{5} + 10 \, b^{8} c d x^{3} - 10 \, a b^{7} d^{2} x^{3} + 15 \, b^{8} c^{2} x - 60 \, a b^{7} c d x + 45 \, a^{2} b^{6} d^{2} x}{15 \, b^{10}} \end{align*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate(x^4*(d*x^2+c)^2/(b*x^2+a)^2,x, algorithm="giac")

[Out]

-1/2*(3*a*b^2*c^2 - 10*a^2*b*c*d + 7*a^3*d^2)*arctan(b*x/sqrt(a*b))/(sqrt(a*b)*b^4) + 1/2*(a*b^2*c^2*x - 2*a^2
*b*c*d*x + a^3*d^2*x)/((b*x^2 + a)*b^4) + 1/15*(3*b^8*d^2*x^5 + 10*b^8*c*d*x^3 - 10*a*b^7*d^2*x^3 + 15*b^8*c^2
*x - 60*a*b^7*c*d*x + 45*a^2*b^6*d^2*x)/b^10